Exploratory analysis of Perceptions about Science and Open Science¶
This analysis revolves around the undergraduate thesis carried out by Franco Sebastián Benítez, under the supervision of Débora Burin and Lucas Cuenya, from the School of Psychology of the University of Buenos Aires.
As set in our preregistration, we are checking the following aspects:
1) Check for exclusion criteria in the demographic data, and in the completion rate.
2) Describe the sample’s demographic characteristics.
3) Analyse the total percentage of “yes” responses to belief in crisis. Analyse as a function of career stage and methodological approach.
4) Qualitative analysis of open field response to belief in crisis.
5) Percentage of agreement with each, and combined, statements about replication crisis, p-value, publication bias. Analyse as a function of career stage and methodological approach.
6) Percentage of agreement with each, and combined, statements about perceived barriers. Analyse as a function of career stage and methodological approach.
7) Percentage of agreement with each, and combined, statements about attitudes against adopting open science practices. Analyse as a function of career stage and methodological approach.
8) Qualitative analysis of open field response to attitudes about barriers against adopting open science practices.
In [179]:
from IPython.core.display import HTML
HTML("""
<style>
.output_png {
display: table-cell;
text-align: center;
vertical-align: middle;
}
</style>
""")
Out[179]:
Loading the necessary libraries¶
In [1]:
import pandas as pd # data wrangling
import matplotlib.pyplot as plt # plotting
import seaborn as sns # plotting 2.0
import re # format text
from wordcloud import WordCloud, STOPWORDS # text analysis
from unidecode import unidecode # remove accents
import numpy as np # numerical manipulations
Loading the dataset¶
In [2]:
df = pd.read_csv("../data/Percepciones sobre ciencia y ciencia abierta.csv")
In [3]:
pd.set_option('max_columns', None) # display all columns
pd.set_option('display.max_colwidth', None) # display full column information
In [4]:
df.head()
Out[4]:
In [5]:
df.shape
Out[5]:
As we see, the data contains 95 rows and 53 columns.
Now, let's rename the columns to make it easier to manipulate and plot.
In [6]:
column_names = {"Timestamp": "timestamp",
"¿Acepta participar?": "consent",
"Edad (años)": "age",
"Nivel educativo alcanzado": "education",
"Área/s de investigación": "area",
"¿Ha participado en un proyecto de investigación (v. g., UBACyT, CONICET) en los últimos 5 años?": "project",
"¿Ha publicado en una revista indexada con referato (v. g., Scopus, Scimago, Scielo) en los últimos 5 años?": "journal",
"Marque su posición actual en la Facultad de Psicología de la UBA": "position",
"¿Qué tipo de metodología suele predominar en sus estudios?": "methodology",
"¿Cree que hay una crisis en la ciencia?": "belief",
"Si su respuesta a la pregunta anterior ha sido “Sí”, señale por qué cree que hay una crisis en la ciencia": "belief_comments",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [Gran parte de los investigadores solamente publica los estudios en donde obtienen diferencias estadísticamente significativas]": "Gran parte de los investigadores solamente publica los estudios en donde obtienen diferencias estadísticamente significativas",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [A las revistas científicas no les interesa publicar resultados que no son estadísticamente significativos]": "A las revistas científicas no les interesa publicar resultados que no son estadísticamente significativos",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [Cuando leo un artículo tengo mayor confianza en la calidad del estudio si los resultados son estadísticamente significativos]": "Cuando leo un artículo tengo mayor confianza en la calidad del estudio si los resultados son estadísticamente significativos",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [Que un resultado no fuera estadísticamente significativo sería un criterio para no seguir investigando esas variables]": "Que un resultado no fuera estadísticamente significativo sería un criterio para no seguir investigando esas variables",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [Una conclusión científica (por ejemplo, si un tratamiento es mejor que otro) debe estar basado en si el p-valor es o no es estadísticamente significativo]": "Una conclusión científica (por ejemplo, si un tratamiento es mejor que otro) debe estar basado en si el p-valor es o no es estadísticamente significativo",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [Un resultado estadísticamente significativo es un resultado importante]": "Un resultado estadísticamente significativo es un resultado importante",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [El valor de p < .05 confirma que el hallazgo será útil para la comunidad científica]": "El valor de p < .05 confirma que el hallazgo será útil para la comunidad científica",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [El valor de p = .001 confirma que el tamaño del efecto ha sido grande]": "El valor de p = .001 confirma que el tamaño del efecto ha sido grande",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [Cuando los resultados de diferentes investigadores son contradictorios entonces sería necesario llevar a cabo un estudio de replicación sobre la misma temática]": "Cuando los resultados de diferentes investigadores son contradictorios entonces sería necesario llevar a cabo un estudio de replicación sobre la misma temática",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [La realización de estudios de replicación es necesaria para el avance de la ciencia]": "La realización de estudios de replicación es necesaria para el avance de la ciencia",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [Si los resultados sobre una determinada temática son unánimes por parte de diferentes equipos de investigación entonces los estudios de replicación no son necesarios]": "Si los resultados sobre una determinada temática son unánimes por parte de diferentes equipos de investigación entonces los estudios de replicación no son necesarios",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [Los estudios de replicación solamente tienen sentido cuando no se detectan diferencias estadísticamente significativas en el estudio original]": "Los estudios de replicación solamente tienen sentido cuando no se detectan diferencias estadísticamente significativas en el estudio original",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [No es necesario replicar un estudio cuando ya se detectaron efectos estadísticamente significativos en el estudio original]": "No es necesario replicar un estudio cuando ya se detectaron efectos estadísticamente significativos en el estudio original",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [El principal objetivo de las revistas científicas es publicar hallazgos novedosos]": "El principal objetivo de las revistas científicas es publicar hallazgos novedosos",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [La ciencia avanza más con estudios que plantean hipótesis novedosas que con estudios de replicación de otras investigaciones]": "La ciencia avanza más con estudios que plantean hipótesis novedosas que con estudios de replicación de otras investigaciones",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [Creo que, en general, los errores que se pueden haber cometido en un estudio científico siempre se detectan y corrigen antes de ser publicados]": "Creo que, en general, los errores que se pueden haber cometido en un estudio científico siempre se detectan y corrigen antes de ser publicados",
"Valore cada una de las siguientes cuestiones relacionadas con su opinión sobre la ciencia [Las tareas de revisión y corrección de los manuscritos que llevan a cabo los revisores de las revistas garantizan de forma fiable la calidad de los resultados científicos]": "Las tareas de revisión y corrección de los manuscritos que llevan a cabo los revisores de las revistas garantizan de forma fiable la calidad de los resultados científicos",
" [Replicar un estudio previo]": "Replicar un estudio previo (1)",
" [Pre-registrar un estudio antes de iniciar la recolección de datos]": "Pre-registrar un estudio antes de iniciar la recolección de datos (1)",
" [Poner datos recolectados a disposición mediante plataformas o repositorios públicos.]": "Poner datos recolectados a disposición mediante plataformas o repositorios públicos (1)",
" [Compartir los métodos analíticos (e.g., scripts de análisis) en plataformas o repositorios públicos]": "Compartir los métodos analíticos (e.g., scripts de análisis) en plataformas o repositorios públicos (1)",
" [Poner los materiales (e.g., cuestionarios, procedimientos) a disposición mediante plataformas o repositorios públicos]": "Poner los materiales (e.g., cuestionarios, procedimientos) a disposición mediante plataformas o repositorios públicos (1)",
" [Reportar toda la información necesaria detalladamente para que otras personas puedan replicar mi estudio]": "Reportar toda la información necesaria detalladamente para que otras personas puedan replicar mi estudio (1)",
" [Compartir una pre-impresión (pre-print) publicándola en un repositorio público o semi-público]": "Compartir una pre-impresión (pre-print) publicándola en un repositorio público o semi-público (1)",
" [Usar revisión por pares abierta]": "Usar revisión por pares abierta (1)",
" [Publicar en revistas de acceso abierto]": "Publicar en revistas de acceso abierto (1)",
" [Adaptar un test psicométrico]": "Adaptar un test psicométrico (1)",
" [Participar en proyectos abiertos y colaborativos a gran escala]": "Participar en proyectos abiertos y colaborativos a gran escala (1)",
" [Replicar un estudio previo].1": "Replicar un estudio previo (2)",
" [Pre-registrar un estudio antes de iniciar la recolección de datos].1": "Pre-registrar un estudio antes de iniciar la recolección de datos (2)",
" [Poner datos recolectados a disposición mediante plataformas o repositorios públicos]": "Poner datos recolectados a disposición mediante plataformas o repositorios públicos (2)",
" [Compartir los métodos analíticos (e.g., scripts de análisis) en plataformas o repositorios públicos].1": "Compartir los métodos analíticos (e.g., scripts de análisis) en plataformas o repositorios públicos (2)",
" [Poner los materiales (e.g., cuestionarios, procedimientos) a disposición mediante plataformas o repositorios públicos].1": "Poner los materiales (e.g., cuestionarios, procedimientos) a disposición mediante plataformas o repositorios públicos (2)",
" [Reportar toda la información necesaria detalladamente para que otras personas puedan replicar mi estudio].1": "Reportar toda la información necesaria detalladamente para que otras personas puedan replicar mi estudio (2)",
" [Compartir una pre-impresión (pre-print) publicándola en un repositorio de confianza]": "Compartir una pre-impresión (pre-print) publicándola en un repositorio de confianza (2)",
" [Usar revisión por partes abierta]": "Usar revisión por partes abierta (2)",
" [Publicar en revistas de acceso abierto].1": "Publicar en revistas de acceso abierto (2)",
" [Adaptar un test psicométrico].1": "Adaptar un test psicométrico (2)",
" [Participar en proyectos abiertos y colaborativos a gran escala].1": "Participar en proyectos abiertos y colaborativos a gran escala (2)",
"¿Cuáles le parecen que son las mayores barreras para la aceptación y puesta en práctica de prácticas de ciencia abierta en su campo y/o lugar de trabajo?": "barriers_1",
"Si aplica, por favor describa brevemente qué barreras ha experimentado para incorporar o mantener prácticas de ciencia abierta": "barriers_2",
"Por último, si posee alguna idea o comentario respecto a esta encuesta o al tema que aborda, por favor escribalo brevemente a continuación ": "comments"
}
df.rename(columns = column_names, inplace=True)
df.timestamp = pd.to_datetime(df.timestamp)
df.age = [re.sub(r"[a-zñ]+", "", i) for i in df["age"]] # remove " años" in values, to convert them in integers
df.age = pd.to_numeric(df.age, downcast="integer")
Let's check all the types.
In [7]:
df.dtypes
Out[7]:
Let's check all the renamings.
In [8]:
df.columns
Out[8]:
Analysis¶
1) Check for exclusion criteria in the demographic data, and in the completion rate¶
To be considered "researchers", our respondents should either have participated in a research project or have publicated in a scientific journal. In both cases, in a maximum range of five years.
In [9]:
exclusion_criteria = df[["project", "journal"]]
excluded = []
for (i, (a, b)) in enumerate(exclusion_criteria.itertuples(index=False), start=0):
if a == "No" and b == "No":
print(f"Participant n° {i} will be excluded")
excluded.append(i)
df = df.drop(excluded, axis=0)
In [10]:
df.shape
Out[10]:
As we can observe, now we have five participants (rows) less.
2) Describe the sample’s demographic characteristics¶
2.1) Education¶
We see the data.
In [11]:
df["education"]
Out[11]:
And we plot them directly, creating a function that we are going to reuse later.
In [12]:
def horizontal_bar_plot(serie, title=None):
y = serie.value_counts(ascending=True)
fig, ax = plt.subplots(figsize=(20,15))
ax.tick_params(labelsize=20)
ax.barh(y.index, y, height=0.75)
ax.set_title(title, fontsize=20)
ax.set_xlabel("Researchers", fontsize=20)
ax.set_ylabel("Frecuency", fontsize=20)
_, xmax = plt.xlim()
plt.xlim(0, xmax)
for i, v in enumerate(y):
ax.text(v + .5, i, str(v), fontsize=20, ha='left', va='center')
plt.show()
In [13]:
horizontal_bar_plot(df["education"], "Education")
Let's go to clean the data a little. We will group each category into four large groups: "Doctorado", "Licenciatura", "Especialización", and "Maestría" regardless of whether it is ongoing or completed. In addition, since we only have one "Postdoctorado", the same will be located in "Doctorado".
First, let's go to separate our different categories by semicolon. For that, we create a new dataframe specific to our variable. The same contains more rows because of the previous rows were expanded by semicolon. The variable, in this case, will be called "education_df".
Note: We add the "belief" variable because we are going to use this variable later.
In [14]:
education_df = df[["education",
"belief",
'Gran parte de los investigadores solamente publica los estudios en donde obtienen diferencias estadísticamente significativas',
'A las revistas científicas no les interesa publicar resultados que no son estadísticamente significativos',
'Cuando leo un artículo tengo mayor confianza en la calidad del estudio si los resultados son estadísticamente significativos',
'Que un resultado no fuera estadísticamente significativo sería un criterio para no seguir investigando esas variables',
'Una conclusión científica (por ejemplo, si un tratamiento es mejor que otro) debe estar basado en si el p-valor es o no es estadísticamente significativo',
'Un resultado estadísticamente significativo es un resultado importante',
'El valor de p < .05 confirma que el hallazgo será útil para la comunidad científica',
'El valor de p = .001 confirma que el tamaño del efecto ha sido grande',
'Cuando los resultados de diferentes investigadores son contradictorios entonces sería necesario llevar a cabo un estudio de replicación sobre la misma temática',
'La realización de estudios de replicación es necesaria para el avance de la ciencia',
'Si los resultados sobre una determinada temática son unánimes por parte de diferentes equipos de investigación entonces los estudios de replicación no son necesarios',
'Los estudios de replicación solamente tienen sentido cuando no se detectan diferencias estadísticamente significativas en el estudio original',
'No es necesario replicar un estudio cuando ya se detectaron efectos estadísticamente significativos en el estudio original',
'El principal objetivo de las revistas científicas es publicar hallazgos novedosos',
'La ciencia avanza más con estudios que plantean hipótesis novedosas que con estudios de replicación de otras investigaciones',
'Creo que, en general, los errores que se pueden haber cometido en un estudio científico siempre se detectan y corrigen antes de ser publicados',
'Las tareas de revisión y corrección de los manuscritos que llevan a cabo los revisores de las revistas garantizan de forma fiable la calidad de los resultados científicos',
'Replicar un estudio previo (1)',
'Pre-registrar un estudio antes de iniciar la recolección de datos (1)',
'Poner datos recolectados a disposición mediante plataformas o repositorios públicos (1)',
'Compartir los métodos analíticos (e.g., scripts de análisis) en plataformas o repositorios públicos (1)',
'Poner los materiales (e.g., cuestionarios, procedimientos) a disposición mediante plataformas o repositorios públicos (1)',
'Reportar toda la información necesaria detalladamente para que otras personas puedan replicar mi estudio (1)',
'Compartir una pre-impresión (pre-print) publicándola en un repositorio público o semi-público (1)',
'Usar revisión por pares abierta (1)',
'Publicar en revistas de acceso abierto (1)',
'Adaptar un test psicométrico (1)',
'Participar en proyectos abiertos y colaborativos a gran escala (1)',
'Replicar un estudio previo (2)',
'Pre-registrar un estudio antes de iniciar la recolección de datos (2)',
'Poner datos recolectados a disposición mediante plataformas o repositorios públicos (2)',
'Compartir los métodos analíticos (e.g., scripts de análisis) en plataformas o repositorios públicos (2)',
'Poner los materiales (e.g., cuestionarios, procedimientos) a disposición mediante plataformas o repositorios públicos (2)',
'Reportar toda la información necesaria detalladamente para que otras personas puedan replicar mi estudio (2)',
'Compartir una pre-impresión (pre-print) publicándola en un repositorio de confianza (2)',
'Usar revisión por partes abierta (2)',
'Publicar en revistas de acceso abierto (2)',
'Adaptar un test psicométrico (2)',
'Participar en proyectos abiertos y colaborativos a gran escala (2)'
]].assign(education=df["education"].str.split(";")).explode("education").reset_index(drop=True)
Then, we have the next new dataframe:
In [15]:
education_df
Out[15]:
Now, let's replace, creating a new function that we are using later.
In [16]:
def replace_columns(df, column_to_replace, original, replaced):
"""
Returns: a df and its columns with strings replaced.
"""
for col in [column_to_replace]:
df[col] = df[col].str.replace(original, replaced)
replace_columns(education_df, "education", "Licenciado", "Licenciatura")
replace_columns(education_df, "education", "Estudiante de licenciatura", "Licenciatura")
replace_columns(education_df, "education", "estudiante de doctorado", "Doctorado")
replace_columns(education_df, "education", "Doctorado en curso", "Doctorado")
replace_columns(education_df, "education", "Estudiante de Doctorado", "Doctorado")
replace_columns(education_df, "education", "En proceso de tesis de maestría", "Maestría")
replace_columns(education_df, "education", "Postdoctorado", "Doctorado")
And we plot, creating a function for later, this time vertically.
In [17]:
def vertical_bar_plot(df, serie, title=None):
plt.figure(figsize=(15,10))
ax = df.groupby(serie).size().sort_values(ascending=False).plot.bar()
ax.set_title(title)
ax.set_xlabel("Researchers")
ax.set_ylabel("Frecuency")
def add_value_labels(ax, spacing=5):
for rect in ax.patches:
y_value = rect.get_height()
x_value = rect.get_x() + rect.get_width() / 2
space = spacing
va = 'bottom'
if y_value < 0:
space *= -1
va = 'top'
label = y_value
ax.annotate(
label,
(x_value, y_value),
xytext=(0, space),
textcoords="offset points",
ha='center',
va=va)
add_value_labels(ax)
plt.xticks(rotation=0)
In [18]:
vertical_bar_plot(education_df, "education", "Education")
Thereby, we have that most respondents have or are getting a doctorate grade, followed by those that have or are getting a licentiate degree.
2.2) Research area¶
First, we create a function that allow us to: 1) remove accents and uppercases, 2) join text from all rows in serie, and 3) drop NA's in case of having.
In [19]:
def extract_text(serie):
serie.dropna(inplace=True)
serie = serie.apply(unidecode).str.casefold().reset_index(drop=True)
text = " ".join(text for text in serie)
return text
Now, we create a function to plot. (We are going to reuse this function later.)
In [20]:
def make_worldcloud(data, stopwords=None):
wordcloud=WordCloud(
background_color='white',
stopwords=stopwords,
max_words=800,
max_font_size=800,
width=800, height=800,
random_state=0
).generate(data)
plt.figure(figsize = (12,15), facecolor=None)
plt.imshow(wordcloud)
plt.axis("off")
plt.tight_layout(pad = 0)
return plt.show()
Finally, we extract text from a serie and plot it.
In [21]:
make_worldcloud(extract_text(df["area"]))
We run again adding a list of stopwords as argument.
In [22]:
make_worldcloud(extract_text(df["area"]), ["psicologia", "y", "la", "de"])
We can see that neuroscience, neuropsychology, social psychology ("social"), developmentental psychology ("del desarrollo"), clinical psychology ("clinica"), and health psychology ("salud"), seem to be the most frequent areas in our sample.
2. 3) Position¶
Let's reuse the function to create a horizontal bar plot.
In [23]:
horizontal_bar_plot(df["position"])
As it was done previously, we store the new values to a new dataframe. In this case, it will be called "position_df".
In [24]:
position_df = df[["position", "belief"]].assign(position=df["position"].str.split(";")).explode("position").reset_index(drop=True)
We plot the data.
In [25]:
horizontal_bar_plot(position_df["position"])
Better, but it still does not look very good.
We group the unique values in a category we will call "Other".
In [26]:
position_df.loc[position_df['position'].isin((position_df['position'].value_counts()[position_df['position'].value_counts() == 1]).index), 'position'] = 'Other'
We plot the data again.
In [27]:
horizontal_bar_plot(position_df["position"])
We can see that most respondents have the rol of "Ayudante de Trabajos Prácticos", by its acronym "ATP".
2.4) Methodology¶
Now, let's see the type of methodology that predominates in our sample of researchers. Directly, we plot the methodology serie from our main dataframe.
In [28]:
df["methodology"].value_counts().plot.pie(figsize = (10,5), autopct='%1.0f%%', fontsize=20)
Out[28]:
As it can be seen, most respondents consider themselves predominantly practising a quantitative approach, being approximately a quarter those that strictly practise a qualitative methodology.
2.5) Age¶
As to respondents' age, first we extract statistics and, next, plot the data directly.
In [29]:
df["age"].describe()
Out[29]:
In [30]:
sns.distplot(df['age'])
Out[30]:
As we observe, most respondents are located in the 20-40 age group.
3) Analyse the total percentage of “yes” responses to belief in crisis. Analyse as a function of career stage and methodological approach¶
3.1) Belief¶
Regarding the belief in crisis variable, we directly plot the data.
In [31]:
df["belief"].value_counts().plot.pie(title="Do you believe that there is a crisis in science?",
figsize = (10,5),
autopct = "%1.0f%%",
fontsize = 20)
Out[31]:
As we can see, the results are almost divided. Other surveys, such as the Baker (2016) survey, perceive higher results as to "yes" percentages.
Below we will be seeing how those respondents who said "yes" justify their answers.
3.2) Belief as a function of career stage¶
Now, we group the belief in crisis by career stage (researchers' education). We plot the data directly, enlarging the plot size.
In [32]:
plt.figure(figsize=(15,10))
sns.countplot(x='education', hue='belief', data=education_df)
Out[32]:
We don't observe a big difference analyzing by career stage.
We can also analyze by position in college.
In [33]:
plt.figure(figsize=(15,10))
sns.countplot(x='position', hue='belief', data=position_df)
Out[33]:
In this case, ATP's seem to believe a little more in crisis in science.
3.3) Belief as a function of methodological approach¶
In [34]:
sns.countplot(x="methodology", hue="belief", data=df)
Out[34]:
Let's take a closer look.
In [35]:
x, y = 'methodology', "belief"
df1 = df.groupby(x)[y].value_counts(normalize=True)
df1 = df1.mul(100)
df1 = df1.rename('percent').reset_index()
g = sns.catplot(x=x,
y='percent',
hue=y,
kind='bar',
data=df1,
height=5,
hue_order=["No", "Sí"],
aspect=2)
for ax in g.axes.ravel():
for p in ax.patches:
ax.annotate(format(p.get_height(), '.2f'), (p.get_x() + p.get_width() / 2.,
p.get_height()), ha = 'center', va = 'center', xytext = (0, 10),
textcoords = 'offset points')
Considering methodology, we can see that there are not substantial differences between usage of mixted and quantitative methodology as to belief in a crisis in science. However, there is a small difference between those that use predominatly qualitative methodology. That is, qualitative researchers seem to believe more in a crisis in science.
Likewise, it is worth noting that the results between those using predominantly quantitative methodology are divided, given that the called "replicability crisis" has a lot to do with statistical problems, such as huge confidence in p-value and null hypothesis testing, and statistical fallacies.
3.4) Belief as a function of age¶
In [36]:
sns.boxplot(data=df, x='belief', y='age')
Out[36]:
In regard to respondents' age, there is not a clear correlation with the belief in crisis in science variable.
4) Qualitative analysis of open field response to belief in crisis¶
Now, let's go to reuse the make_worldcloud() and extract_text() functions, built previously, to analyse justifications to "yes" answers to belief in crisis.
In [37]:
make_worldcloud(extract_text(df["belief_comments"]))
It doesn't look very good.
We add stopwords.
In [38]:
make_worldcloud(extract_text(df["belief_comments"]), ["que", "la", "y", "por", "porque", "en", "el", "para", "de", "se", "lo"
, "a", "como", "tambien", "los", "las", "del", "hay", "una", "pero",
"esta", "ciencia", "investigadores", "investigacion", "no", "su", "un",
"con", "al", "e", "mucha", "crisis", "o", "es", "cientifico"])
As now we can see, having filtered vague words, there are some words more frequent among those that answered "yes" to the question "Do you believe there are a crisis in science?". Such words more frequent are replicability ("replicabilidad"), system ("sistema"), absence ("falta"), and quality ("calidad").
Now let's see how many comments mentioned the "replicability" word as a cause of the crisis. First we create a new dataframe removing accents and uppercase of the comments.
In [39]:
belief_comments_df = df["belief_comments"].dropna().apply(unidecode).str.casefold().to_frame().reset_index(drop=True)
In [40]:
belief_comments_df["belief_comments"].count()
Out[40]:
In [41]:
l = [text for text in belief_comments_df["belief_comments"]]
In [42]:
for (i, text) in enumerate(l):
print('Comment %d:\n'%(i+1), text)
We add the comments mentioning "replica" to a list. We use "replica" to catch both "replicacion" and "replicabilidad".
In [43]:
l = [text for text in belief_comments_df["belief_comments"] if "replica" in text]
We print the data formatted.
In [44]:
for (i, text) in enumerate(l):
print('Comment %d:\n'%(i+1), text)
As we can see, only seven comments mention "replicabilidad" or "replicación" as causes of crisis in science.
5) Percentage of agreement with each, and combined, statements about replication crisis, p-value, publication bias. Analyse as a function of career stage and methodological approach¶
In this point, let's go to create tidy tables that allow us to plot easily. These data, given that are provided for a likert scale, will be plotted with the HH package in R. The plotting code will be available in the R Script.
First, let's create a function that builds a table from a column provided.
In [45]:
def make_table(column):
"""
Returns: a table.
"""
df = column.value_counts().to_frame().reset_index()
df = df.T
df.columns = df.iloc[0]
df = df.drop(df.index[0])
return df
Second, we apply the previous function to create tables iteratively from an index to other.
In [46]:
def tables_by_columns(index_col_1, index_col_2):
"""
Returns: a list containing different tables created from indexes.
"""
l = []
for (i, column) in enumerate(df):
if (i >= index_col_1 and i <= index_col_2):
l.append(make_table(df[column]))
return l
Third, we combinate the tables creating only one.
In [47]:
def combinate_tables(l):
"""
Returns: a new dataframe.
"""
df = pd.concat(l)
df.reset_index(level=0, inplace=True)
df = df.rename(columns={'index': 'Measures'})
return df
Finally, we convert numerical values to int, given that by default they are strings.
In [48]:
def convert_to_int(df):
"""
Returns: a df with numerical values converted to int.
"""
for (i, column) in enumerate(df):
if (i > 0):
df[[column]] = df[[column]].astype('Int64').fillna(0)
return df
We are creating three variables relative to:
1) Answers to "Value each one of the following issues concerning your opinion about science". This variable will be called "science".
2) Answers to "Mark the option that best represents your knowledge and experience with each practise in the last five years". This variable will be named "experience".
3) Answers to "Choose the option that best represents how much importance you consider that has each one of the following practises to improve quality and efficiency of research in your research area". This variable will be named "efficiency".
In [49]:
science = convert_to_int(combinate_tables(tables_by_columns(11, 27)))
experience = convert_to_int(combinate_tables(tables_by_columns(28, 38)))
efficiency = convert_to_int(combinate_tables(tables_by_columns(39, 49)))
Let's see the results.
In [50]:
science
Out[50]:
We sort labels.
In [51]:
science = science[["Measures", "Nada de acuerdo", "Algo de acuerdo", "No sé", "Bastante de acuerdo", "Muy de acuerdo"]]
experience = experience[["Measures",
"Nunca lo hice, ni escuché hablar de eso",
"Nunca lo hice, pero escuché hablar de eso",
"Nunca lo hice, pero sé cómo hacerlo",
"Lo intenté, pero no lo completé",
"Sí, lo estoy haciendo",
"Sí, lo he hecho"]]
efficiency = efficiency[["Measures",
"Nada importante",
"Poco importante",
"No sé",
"Medianamente importante",
"Muy importante"]]
In [52]:
science
Out[52]:
In [53]:
experience
Out[53]:
In [54]:
efficiency
Out[54]:
We export (to be plotted from R and the HH package).
In [55]:
science.to_csv(r'../data/cleaned/science.csv', index=False)
experience.to_csv(r'../data/cleaned/experience.csv', index=False)
efficiency.to_csv(r'../data/cleaned/efficiency.csv', index=False)
And we import the resulting plots.
(1) Science:
In [56]:
#
(2) Experience:
In [57]:
#
(3) Efficiency:
In [58]:
#
The plots in the best quality can be found in https://github.com/francosbenitez/thesis/tree/master/images.
In [59]:
#Now, we create and export tables grouped by career stage and methodological approach.
#First, we create a function that receives a column and a group and transform them in a dataframe table.
In [60]:
# def make_table_by_group(column, group):
# """
# Returns: a table by group in dataframe format.
# """
# new_df = df.groupby(column)[group].value_counts().unstack().fillna(0).reset_index()
# new_df = new_df.T
# new_df.columns = new_df.iloc[0]
# new_df = new_df.drop(new_df.index[0])
# new_df.columns.name = None
# new_df = new_df.reset_index()
# return new_df
In [61]:
#Then, we create iteratively tables by index, but, unlike before, this time will be by group.
In [62]:
# def tables_by_columns_and_group(index_col_1, index_col_2, group):
# """
# Returns: a dictionary with a variable name as key and a table as value.
# """
# d = {}
# column_names = {}
# for (i, column) in enumerate(df):
# if (i >= index_col_1 and i <= index_col_2):
# column_names[column] = make_table_by_group(df[column], group)
# d.update(column_names)
# return d
In [63]:
#We create a function to export iteratively.
In [64]:
# def export_table(table, name):
# """
# Returns: different tables in csv format.
# """
# for (i, (k, v)) in enumerate(table.items()):
# v.to_csv("../data/cleaned/"+name+"_"+str(i)+".csv")
In [65]:
#We export.
In [66]:
# export_table(tables_by_columns_and_group(11, 27, "education"), "science_by_education")
# export_table(tables_by_columns_and_group(28, 38, "education"), "experience_by_education")
# export_table(tables_by_columns_and_group(39, 49, "education"), "efficiency_by_education")
# export_table(tables_by_columns_and_group(11, 27, "methodology"), "science_by_methodology")
# export_table(tables_by_columns_and_group(28, 38, "methodology"), "experience_by_methodology")
# export_table(tables_by_columns_and_group(39, 49, "methodology"), "efficiency_by_methodology")
In [67]:
def science_by_methodology(serie):
x,y = 'methodology', serie
df1 = df.groupby(x)[y].value_counts(normalize=True)
df1 = df1.mul(100)
df1 = df1.rename('percent').reset_index()
g = sns.catplot(x=x,
y='percent',
hue=y,
kind='bar',
data=df1,
height=5,
aspect=2,
hue_order = ["Nada de acuerdo", "Algo de acuerdo", "No sé", "Bastante de acuerdo", "Muy de acuerdo"])
g.ax.set_ylim(0,100)
for ax in g.axes.ravel():
for p in ax.patches:
ax.annotate(format(p.get_height(), '.2f'), (p.get_x() + p.get_width() / 2.,
p.get_height()), ha = 'center', va = 'center', xytext = (0, 10),
textcoords = 'offset points')
g._legend.remove()
plt.legend(loc='upper right')
plt.title(serie + "\n" + "\n")
In [68]:
science_by_methodology("Cuando los resultados de diferentes investigadores son contradictorios entonces sería necesario llevar a cabo un estudio de replicación sobre la misma temática")
In [69]:
science_by_methodology("Cuando los resultados de diferentes investigadores son contradictorios entonces sería necesario llevar a cabo un estudio de replicación sobre la misma temática")
In [70]:
science_by_methodology("La realización de estudios de replicación es necesaria para el avance de la ciencia")
In [71]:
science_by_methodology("A las revistas científicas no les interesa publicar resultados que no son estadísticamente significativos")
In [72]:
science_by_methodology("Gran parte de los investigadores solamente publica los estudios en donde obtienen diferencias estadísticamente significativas")
In [73]:
science_by_methodology("El principal objetivo de las revistas científicas es publicar hallazgos novedosos")
In [74]:
science_by_methodology("Un resultado estadísticamente significativo es un resultado importante")
In [75]:
science_by_methodology("Una conclusión científica (por ejemplo, si un tratamiento es mejor que otro) debe estar basado en si el p-valor es o no es estadísticamente significativo")
In [76]:
science_by_methodology("Cuando leo un artículo tengo mayor confianza en la calidad del estudio si los resultados son estadísticamente significativos")
In [77]:
science_by_methodology("El valor de p = .001 confirma que el tamaño del efecto ha sido grande")
In [78]:
science_by_methodology("Las tareas de revisión y corrección de los manuscritos que llevan a cabo los revisores de las revistas garantizan de forma fiable la calidad de los resultados científicos")
In [79]:
science_by_methodology("El valor de p < .05 confirma que el hallazgo será útil para la comunidad científica")
In [80]:
science_by_methodology("La ciencia avanza más con estudios que plantean hipótesis novedosas que con estudios de replicación de otras investigaciones")
In [81]:
science_by_methodology("Si los resultados sobre una determinada temática son unánimes por parte de diferentes equipos de investigación entonces los estudios de replicación no son necesarios")
In [82]:
science_by_methodology("Creo que, en general, los errores que se pueden haber cometido en un estudio científico siempre se detectan y corrigen antes de ser publicados")
In [83]:
science_by_methodology("Los estudios de replicación solamente tienen sentido cuando no se detectan diferencias estadísticamente significativas en el estudio original")
In [84]:
science_by_methodology("No es necesario replicar un estudio cuando ya se detectaron efectos estadísticamente significativos en el estudio original")
In [85]:
science_by_methodology("Que un resultado no fuera estadísticamente significativo sería un criterio para no seguir investigando esas variables")
In [86]:
def science_by_education(serie):
x,y = 'education', serie
df1 = education_df.groupby(x)[y].value_counts(normalize=True)
df1 = df1.mul(100)
df1 = df1.rename('percent').reset_index()
g = sns.catplot(x=x,
y='percent',
hue=y,
kind='bar',
data=df1,
#orient="h",
height=5,
aspect=2,
hue_order = ["Nada de acuerdo", "Algo de acuerdo", "No sé", "Bastante de acuerdo", "Muy de acuerdo"])
g.ax.set_ylim(0,100)
for ax in g.axes.ravel():
for p in ax.patches:
ax.annotate(format(p.get_height(), '.2f'), (p.get_x() + p.get_width() / 2.,
p.get_height()), ha = 'center', va = 'center', xytext = (0, 10),
textcoords = 'offset points')
g._legend.remove()
plt.legend(loc='upper right')
plt.title(serie + "\n" + "\n")
In [87]:
science_by_education("Cuando los resultados de diferentes investigadores son contradictorios entonces sería necesario llevar a cabo un estudio de replicación sobre la misma temática")
In [88]:
science_by_education("Cuando los resultados de diferentes investigadores son contradictorios entonces sería necesario llevar a cabo un estudio de replicación sobre la misma temática")
In [89]:
science_by_education("La realización de estudios de replicación es necesaria para el avance de la ciencia")
In [90]:
science_by_education("A las revistas científicas no les interesa publicar resultados que no son estadísticamente significativos")
In [91]:
science_by_education("Gran parte de los investigadores solamente publica los estudios en donde obtienen diferencias estadísticamente significativas")
In [92]:
science_by_education("El principal objetivo de las revistas científicas es publicar hallazgos novedosos")
In [93]:
science_by_education("Un resultado estadísticamente significativo es un resultado importante")
In [94]:
science_by_education("Una conclusión científica (por ejemplo, si un tratamiento es mejor que otro) debe estar basado en si el p-valor es o no es estadísticamente significativo")
In [95]:
science_by_education("Cuando leo un artículo tengo mayor confianza en la calidad del estudio si los resultados son estadísticamente significativos")
In [96]:
science_by_education("El valor de p = .001 confirma que el tamaño del efecto ha sido grande")
In [97]:
science_by_education("Las tareas de revisión y corrección de los manuscritos que llevan a cabo los revisores de las revistas garantizan de forma fiable la calidad de los resultados científicos")
In [98]:
science_by_education("El valor de p < .05 confirma que el hallazgo será útil para la comunidad científica")
In [99]:
science_by_education("La ciencia avanza más con estudios que plantean hipótesis novedosas que con estudios de replicación de otras investigaciones")
In [100]:
science_by_education("Si los resultados sobre una determinada temática son unánimes por parte de diferentes equipos de investigación entonces los estudios de replicación no son necesarios")
In [101]:
science_by_education("Creo que, en general, los errores que se pueden haber cometido en un estudio científico siempre se detectan y corrigen antes de ser publicados")
In [102]:
science_by_education("Los estudios de replicación solamente tienen sentido cuando no se detectan diferencias estadísticamente significativas en el estudio original")
In [103]:
science_by_education("No es necesario replicar un estudio cuando ya se detectaron efectos estadísticamente significativos en el estudio original")
In [104]:
science_by_education("Que un resultado no fuera estadísticamente significativo sería un criterio para no seguir investigando esas variables")
In [105]:
def experience_by_methodology(serie):
x,y = 'methodology', serie
df1 = df.groupby(x)[y].value_counts(normalize=True)
df1 = df1.mul(100)
df1 = df1.rename('percent').reset_index()
g = sns.catplot(x=x,
y='percent',
hue=y,
kind='bar',
data=df1,
height=5,
aspect=2,
hue_order = ["Nunca lo hice, ni escuché hablar de eso",
"Nunca lo hice, pero escuché hablar de eso",
"Nunca lo hice, pero sé cómo hacerlo",
"Lo intenté, pero no lo completé",
"Sí, lo estoy haciendo",
"Sí, lo he hecho"])
g.ax.set_ylim(0,100)
for ax in g.axes.ravel():
for p in ax.patches:
ax.annotate(format(p.get_height(), '.2f'), (p.get_x() + p.get_width() / 2.,
p.get_height()), ha = 'center', va = 'center', xytext = (0, 10),
textcoords = 'offset points')
g._legend.remove()
plt.legend(loc='upper right')
plt.title(serie + "\n" + "\n")
In [106]:
experience_by_methodology('Replicar un estudio previo (1)')
In [107]:
experience_by_methodology('Pre-registrar un estudio antes de iniciar la recolección de datos (1)')
In [108]:
experience_by_methodology('Poner datos recolectados a disposición mediante plataformas o repositorios públicos (1)')
In [109]:
experience_by_methodology('Compartir los métodos analíticos (e.g., scripts de análisis) en plataformas o repositorios públicos (1)')
In [110]:
experience_by_methodology('Poner los materiales (e.g., cuestionarios, procedimientos) a disposición mediante plataformas o repositorios públicos (1)')
In [111]:
experience_by_methodology('Reportar toda la información necesaria detalladamente para que otras personas puedan replicar mi estudio (1)')
In [112]:
experience_by_methodology('Compartir una pre-impresión (pre-print) publicándola en un repositorio público o semi-público (1)')
In [113]:
experience_by_methodology('Usar revisión por pares abierta (1)')
In [114]:
experience_by_methodology('Publicar en revistas de acceso abierto (1)')
In [115]:
experience_by_methodology('Adaptar un test psicométrico (1)')
In [116]:
experience_by_methodology('Participar en proyectos abiertos y colaborativos a gran escala (1)')
In [117]:
education_df
Out[117]:
In [118]:
def experience_by_education(serie):
x,y = 'education', serie
df1 = education_df.groupby(x)[y].value_counts(normalize=True)
df1 = df1.mul(100)
df1 = df1.rename('percent').reset_index()
g = sns.catplot(x=x,
y='percent',
hue=y,
kind='bar',
data=df1,
#orient="h",
height=5,
aspect=2,
hue_order = ["Nunca lo hice, ni escuché hablar de eso",
"Nunca lo hice, pero escuché hablar de eso",
"Nunca lo hice, pero sé cómo hacerlo",
"Lo intenté, pero no lo completé",
"Sí, lo estoy haciendo",
"Sí, lo he hecho"])
g.ax.set_ylim(0,100)
for ax in g.axes.ravel():
for p in ax.patches:
ax.annotate(format(p.get_height(), '.2f'), (p.get_x() + p.get_width() / 2.,
p.get_height()), ha = 'center', va = 'center', xytext = (0, 10),
textcoords = 'offset points')
g._legend.remove()
plt.legend(loc='upper right')
plt.title(serie + "\n" + "\n")
In [119]:
experience_by_education('Replicar un estudio previo (1)')
In [120]:
experience_by_education('Pre-registrar un estudio antes de iniciar la recolección de datos (1)')
In [121]:
experience_by_education('Poner datos recolectados a disposición mediante plataformas o repositorios públicos (1)')
In [122]:
experience_by_education('Compartir los métodos analíticos (e.g., scripts de análisis) en plataformas o repositorios públicos (1)')
In [123]:
experience_by_education('Poner los materiales (e.g., cuestionarios, procedimientos) a disposición mediante plataformas o repositorios públicos (1)')
In [124]:
experience_by_education('Reportar toda la información necesaria detalladamente para que otras personas puedan replicar mi estudio (1)')
In [125]:
experience_by_education('Compartir una pre-impresión (pre-print) publicándola en un repositorio público o semi-público (1)')
In [126]:
experience_by_education('Usar revisión por pares abierta (1)')
In [127]:
experience_by_education('Publicar en revistas de acceso abierto (1)')
In [128]:
experience_by_education('Adaptar un test psicométrico (1)')
In [129]:
experience_by_education('Participar en proyectos abiertos y colaborativos a gran escala (1)')
In [130]:
def efficiency_by_methodology(serie):
x,y = 'methodology', serie
df1 = df.groupby(x)[y].value_counts(normalize=True)
df1 = df1.mul(100)
df1 = df1.rename('percent').reset_index()
g = sns.catplot(x=x,
y='percent',
hue=y,
kind='bar',
data=df1,
height=5,
aspect=2,
hue_order = ["Nada importante",
"Poco importante",
"No sé",
"Medianamente importante",
"Muy importante"])
g.ax.set_ylim(0,100)
for ax in g.axes.ravel():
for p in ax.patches:
ax.annotate(format(p.get_height(), '.2f'), (p.get_x() + p.get_width() / 2.,
p.get_height()), ha = 'center', va = 'center', xytext = (0, 10),
textcoords = 'offset points')
g._legend.remove()
plt.legend(loc='upper right')
plt.title(serie + "\n" + "\n")
In [131]:
efficiency_by_methodology('Replicar un estudio previo (2)')
In [132]:
efficiency_by_methodology('Pre-registrar un estudio antes de iniciar la recolección de datos (2)')
In [133]:
efficiency_by_methodology('Poner datos recolectados a disposición mediante plataformas o repositorios públicos (2)')
In [134]:
efficiency_by_methodology('Compartir los métodos analíticos (e.g., scripts de análisis) en plataformas o repositorios públicos (2)')
In [135]:
efficiency_by_methodology('Poner los materiales (e.g., cuestionarios, procedimientos) a disposición mediante plataformas o repositorios públicos (2)')
In [136]:
efficiency_by_methodology('Reportar toda la información necesaria detalladamente para que otras personas puedan replicar mi estudio (2)')
In [137]:
efficiency_by_methodology('Compartir una pre-impresión (pre-print) publicándola en un repositorio de confianza (2)')
In [138]:
efficiency_by_methodology('Usar revisión por partes abierta (2)')
In [139]:
efficiency_by_methodology('Publicar en revistas de acceso abierto (2)')
In [140]:
efficiency_by_methodology('Adaptar un test psicométrico (2)')
In [141]:
efficiency_by_methodology('Participar en proyectos abiertos y colaborativos a gran escala (2)')
In [142]:
def efficiency_by_education(serie):
x,y = 'education', serie
df1 = education_df.groupby(x)[y].value_counts(normalize=True)
df1 = df1.mul(100)
df1 = df1.rename('percent').reset_index()
g = sns.catplot(x=x,
y='percent',
hue=y,
kind='bar',
data=df1,
#orient="h",
height=5,
aspect=2,
hue_order = ["Nada importante",
"Poco importante",
"No sé",
"Medianamente importante",
"Muy importante"])
g.ax.set_ylim(0,100)
for ax in g.axes.ravel():
for p in ax.patches:
ax.annotate(format(p.get_height(), '.2f'), (p.get_x() + p.get_width() / 2.,
p.get_height()), ha = 'center', va = 'center', xytext = (0, 10),
textcoords = 'offset points')
g._legend.remove()
plt.legend(loc='upper right')
plt.title(serie + "\n" + "\n")
In [143]:
efficiency_by_education('Replicar un estudio previo (2)')
In [144]:
efficiency_by_education('Pre-registrar un estudio antes de iniciar la recolección de datos (2)')
In [145]:
efficiency_by_education('Poner datos recolectados a disposición mediante plataformas o repositorios públicos (2)')
In [146]:
efficiency_by_education('Compartir los métodos analíticos (e.g., scripts de análisis) en plataformas o repositorios públicos (2)')
In [147]:
efficiency_by_education('Poner los materiales (e.g., cuestionarios, procedimientos) a disposición mediante plataformas o repositorios públicos (2)')
In [148]:
efficiency_by_education('Reportar toda la información necesaria detalladamente para que otras personas puedan replicar mi estudio (2)')
In [149]:
efficiency_by_education('Compartir una pre-impresión (pre-print) publicándola en un repositorio de confianza (2)')
In [150]:
efficiency_by_education('Usar revisión por partes abierta (2)')
In [151]:
efficiency_by_education('Publicar en revistas de acceso abierto (2)')
In [152]:
efficiency_by_education('Adaptar un test psicométrico (2)')
In [153]:
efficiency_by_education('Participar en proyectos abiertos y colaborativos a gran escala (2)')
6) Percentage of agreement with each, and combined, statements about perceived barriers. Analyse as a function of career stage and methodological approach¶
In [154]:
barriers_1_by_methodology = df[["barriers_1", "methodology"]].assign(barriers_1=df[["barriers_1"]].barriers_1.str.split(";")).explode("barriers_1").reset_index(drop=True)
In [155]:
barriers_1_by_education = df[["barriers_1", "education"]].assign(education=df["education"].str.split(";")).explode("education").reset_index(drop=True)
In [156]:
barriers_1_by_education = barriers_1_by_education[["barriers_1", "education"]].assign(barriers_1=barriers_1_by_education[["barriers_1"]].barriers_1.str.split(";")).explode("barriers_1").reset_index(drop=True)
In [157]:
replace_columns(barriers_1_by_education, "education", "Licenciado", "Licenciatura")
replace_columns(barriers_1_by_education, "education", "Estudiante de licenciatura", "Licenciatura")
replace_columns(barriers_1_by_education, "education", "estudiante de doctorado", "Doctorado")
replace_columns(barriers_1_by_education, "education", "Doctorado en curso", "Doctorado")
replace_columns(barriers_1_by_education, "education", "Estudiante de Doctorado", "Doctorado")
replace_columns(barriers_1_by_education, "education", "En proceso de tesis de maestría", "Maestría")
replace_columns(barriers_1_by_education, "education", "Postdoctorado", "Doctorado")
replace_columns(barriers_1_by_education, "barriers_1", "[0-9.]", "")
In [158]:
g = sns.catplot(y='barriers_1',
col='education',
data=barriers_1_by_education,
kind="count",
aspect=.8,
order=barriers_1_by_education["barriers_1"].value_counts().index)
for ax in g.axes.ravel():
for p in ax.patches:
percentage = '{0:g}'.format(p.get_width())
x = p.get_x() + p.get_width() + 0.02
y = p.get_y() + p.get_height()/2
ax.annotate(percentage, (x, y))
In [159]:
replace_columns(barriers_1_by_methodology, "barriers_1", "[0-9.]", "")
In [160]:
g = sns.catplot(y='barriers_1',
col='methodology',
data=barriers_1_by_methodology,
kind="count",
order=barriers_1_by_methodology["barriers_1"].value_counts().index)
for ax in g.axes.ravel():
for p in ax.patches:
percentage = '{0:g}'.format(p.get_width())
x = p.get_x() + p.get_width() + 0.02
y = p.get_y() + p.get_height()/2
ax.annotate(percentage, (x, y))
As we have already done previously, we create a new dataframe with rows extended by semi-colons.
In [161]:
barriers_1_df = df[["barriers_1"]].assign(barriers_1=df[["barriers_1"]].barriers_1.str.split(";")).explode("barriers_1").reset_index(drop=True)
We replace the numbers at the beginning of each row. (We reuse the replace_columns() function defined previously.)
In [162]:
replace_columns(barriers_1_df, "barriers_1", "[0-9.]", "")
We plot.
In [163]:
horizontal_bar_plot(barriers_1_df["barriers_1"], "Major barriers against adopting open science practices")
Now we plot it adding percents.
In [164]:
def horizontal_bar_plot_with_percents(serie, title=None):
y = serie.value_counts(ascending=True)
fig, ax = plt.subplots(figsize=(20,15))
ax.tick_params(labelsize=20)
ax.barh(y.index, y, height=0.75, color="black")
ax.set_title(title, fontsize=20)
#ax.set_xlabel("Researchers", fontsize=20)
#ax.set_ylabel("Frecuency", fontsize=20)
_, xmax = plt.xlim()
plt.xlim(0, xmax)
plt.xticks([])
for i, v in enumerate(y):
#ax.text(v + .5, i, str('{}%'.format(round(v * 100 / 90), 2)), fontsize=20, ha='left', va='center')
ax.text(v + .5, i, str(format(v * 100 / 90, '.2f')), fontsize=20, ha='left', va='center')
plt.show()
In [165]:
horizontal_bar_plot_with_percents(barriers_1_df["barriers_1"])
In [166]:
### 7) Percentage of agreement with each, and combined, statements about attitudes against adopting open science practices. Analyse as a function of career stage and methodological approach
8) Qualitative analysis of open field response to attitudes about barriers against adopting open science practices¶
In [167]:
make_worldcloud(extract_text(df["barriers_2"]), ["de", "en", "la", "que", "y", "lo", "las", "para", "el", "no",
"por", "a", "e", "los", "mi", "es", "se", "con", "cual", "veces", "muchas",
"mayor", "un"])
In [168]:
barriers_2_df =df["barriers_2"].dropna().apply(unidecode).str.casefold().to_frame().reset_index(drop=True)
In [169]:
barriers_2_df["barriers_2"].count()
Out[169]:
In [170]:
l = [text for text in barriers_2_df["barriers_2"]]
In [171]:
for (i, text) in enumerate(l):
print('Comment %d:\n'%(i+1), text)
9) Other analysis¶
In the end, let's see all comments to the "Finally, if you have any ideas or comments regarding this survey or the topic it covers, please write them briefly below" item.
In [172]:
l = [text for text in df["comments"].dropna()]
In [173]:
for (i, text) in enumerate(l):
print('Comment %d:\n'%(i+1), text)
We add the mean and standart desviation of the tables.
In [174]:
science["M"] = round(science.apply(lambda row: (row["Nada de acuerdo"] * 1) +
(row["Algo de acuerdo"] * 2) +
(row["No sé"] * 3) +
(row["Bastante de acuerdo"] * 4 +
row["Muy de acuerdo"] * 5)
, axis=1) / 90, 2)
efficiency["M"] = round(efficiency.apply(lambda row: (row["Nada importante"] * 1) +
(row["Poco importante"] * 2) +
(row["No sé"] * 3) +
(row["Medianamente importante"] * 4 +
row["Muy importante"] * 5)
, axis=1) / 90, 2)
experience["M"] = round(experience.apply(lambda row: (row["Nunca lo hice, ni escuché hablar de eso"] * 1) +
(row["Nunca lo hice, pero escuché hablar de eso"] * 2) +
(row["Nunca lo hice, pero sé cómo hacerlo"] * 3) +
(row["Lo intenté, pero no lo completé"] * 4 +
row["Sí, lo estoy haciendo"] * 5 +
row["Sí, lo he hecho"] * 6)
, axis=1) / 90, 2)
science["SD"] = science.apply(lambda row: (row["Nada de acuerdo"] * 1) +
(row["Algo de acuerdo"] * 4) +
(row["No sé"] * 9) +
(row["Bastante de acuerdo"] * 16 +
row["Muy de acuerdo"] * 25)
, axis=1) / 90
efficiency["SD"] = efficiency.apply(lambda row: (row["Nada importante"] * 1) +
(row["Poco importante"] * 4) +
(row["No sé"] * 9) +
(row["Medianamente importante"] * 16 +
row["Muy importante"] * 25)
, axis=1) / 90
experience["SD"] = experience.apply(lambda row: (row["Nunca lo hice, ni escuché hablar de eso"] * 1) +
(row["Nunca lo hice, pero escuché hablar de eso"] * 4) +
(row["Nunca lo hice, pero sé cómo hacerlo"] * 9) +
(row["Lo intenté, pero no lo completé"] * 16 +
row["Sí, lo estoy haciendo"] * 25 +
row["Sí, lo he hecho"] * 36)
, axis=1) / 90
science["SD"] = round(np.sqrt(science["SD"] - science["M"]), 2)
experience["SD"] = round(np.sqrt(experience["SD"] - experience["M"]), 2)
efficiency["SD"] = round(np.sqrt(efficiency["SD"] - efficiency["M"]), 2)
In [175]:
science
Out[175]:
In [176]:
experience
Out[176]:
In [177]:
efficiency
Out[177]:
We transform the count values into percentages.
In [178]:
science_cols = ["Nada de acuerdo", "Algo de acuerdo", "No sé", "Bastante de acuerdo", "Muy de acuerdo"]
experience_cols = ["Nunca lo hice, ni escuché hablar de eso",
"Nunca lo hice, pero escuché hablar de eso",
"Nunca lo hice, pero sé cómo hacerlo",
"Lo intenté, pero no lo completé",
"Sí, lo estoy haciendo",
"Sí, lo he hecho"]
efficiency_cols = ["Nada importante",
"Poco importante",
"No sé",
"Medianamente importante",
"Muy importante"]
science[science_cols] = round(science[science_cols].div(science[science_cols].sum(axis=1), axis=0).multiply(100), 1)
experience[experience_cols] = round(experience[experience_cols].div(experience[experience_cols].sum(axis=1), axis=0).multiply(100), 1)
efficiency[efficiency_cols] = round(efficiency[efficiency_cols].div(efficiency[efficiency_cols].sum(axis=1), axis=0).multiply(100), 1)
In [ ]:
science
In [ ]:
experience
In [ ]:
efficiency